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1. Introduction 

The method of large N continuum reduction |2| for SU{N) gauge theory allows for 
the calculation of the infinite volume, infinite N limit of certain physical quantities using 
volumes reduced to a small physical size. Numerical estimates [|l|, || of the physical critical 
size above which continuum reduction holds indicate that this method can be used to 
produce practical results. The chiral condensate |p and the pion decay constant Q| were 
calculated in the large N limit in four dimensions using continuum reduction. In this 
paper, we show that the method can be extended beyond bulk quantities and that it also 
produces reliable results for quantities with space-time dependence such as the heavy quark 
potential, from which the string tension can be extracted. 

A precise calculation of the string tension in three dimensional SU{N) gauge theories 
has been performed with up to 8 on large lattices |^. In this paper, we present a 
complementary calculation with = 47 on 5^ lattices using continuum reduction. The 
calculation of Ref. |^] used correlation functions of smeared Polyakov loops to extract the 
string tension. After extrapolating to A^ = oo and to the continuum, the result was 

= 0.1975 ± 0.0002 - 0.0005 (1.1) 

where g is the gauge coupling. This has to be compared with the analytical calculation 
in Q, namely, ~ 0.1995. Although the two results are not in perfect agreement, the 
main observation is that the approximations used in the analytical calculation are very 
well motivated. 
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The string tension for SU{N) as per the analytical calculation |^ is 



a = = \ q N\ — 

27r 2 J 8^ 



1 



(1.2) 



CA and cp are the quadratic Casimirs in the adjoint and fundamental representation, re- 



2 

spectively. The mass parameter naturally enters the analytical calculation, and the 



^2, 



second factor arises from the Wilson loop operator in the fundamental representation. 

The numerical computation in shows that the agreement gets better as one gets 
closer to A'^ = oo. Like the analytical result, the numerical result also shows a correction 
from N = CO that goes like but the coefficient of is not the same for the numerical 
and the analytical computations. 

Since the numerical and the analytical results are close to each other even for N = 2 
(less than 4%), we can use the analytical formula to get a feel for the finite corrections 
to the infinite A'^ result. Therefore, we expect the finite corrections to be smaller than 



the error in Eqn(l.l) if A^ > 32. 

In this paper, we use continuum reduction ||l|, |2| to directly compute the N = oo limit 
of the string tension by working at large enough A^ so that the finite A^ corrections are 
smaller than the numerical errors. We find that 

= 0.1964 ± 0.0009 (1.3) 

This result and that of ( |1.1| ) are consistent at the level of their one sigma errors. This level 
of agreement is, in turn, consistent with neither the large A^ extrapolation of Ref. nor 
the volume reduction of the present calculation having unexpected errors. While both of 
the numerical results lie below the analytical estimate, the discrepancy is relatively small. 
Thus the numerical evidence that the analytical result is an excellent first approximation 
that captures much of the physics remains strong. 

The paper is organized as follows. We explain how we use smeared Wilson loops to 
compute the string tension in Section ^. The lattice results for the string tension along 
with the continuum extrapolation are also presented in this section. An intermediate step 
in our calculation is the dimensionless ground state string energy m{k). In Section ^, we 
show results for m{k) at one fixed lattice coupling to illustrate its behavior as a function 
of k and how it is used to extract the string tension. We also show that m{k) is unaffected 
by the smearing parameter. We illustrate the extraction of m{k) at one fixed coupling in 
Section ^. Here we show how the smearing parameter affects the overlap with the ground 
state. The main result in this paper is obtained using A^ = 47. We show that the finite A^ 
and finite volume corrections are small at this value of A^ in Section ||. We explain why 
this method is preferred over the Creutz ratio in Section 

2. String tension using Wilson loops and continuum reduction 

Consider SU{N) Yang-Mills theory on a periodic lattice with the standard Wilson gauge 
action. The method of fi] is to measure the string tension using correlations of Polyakov 
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loops with separation t that wind around a space direction. Continuum reduction [||, |2| 
imphes that the large N Yang-Mills theory in a continuum box of size is independent of 
/ as long as / > /c = 1/Tc with Tc being the deconfining temperature. One should be able 
to compute expectation values of Wilson loops of arbitrary size on an l"^ continuum box 
using folded Wilson loops and extract the string tension. To implement this approach to 
the three-dimensional Yang-Mills theory string tension, we use the following procedure: 

• We fix the lattice size to L^. We use L = 5 for the most part and only use L = 4 to 
verify reduction. 

• We fix N so that finite corrections are small. We set N = 47 and show using one 
instance that finite corrections are small at = 47. 

• We pick an appropriate range of lattice coupling b = . 

— b cannot be too small since we have to be away from the bulk transition on the 
lattice associated with the development of gap in the eigenvalue distribution of 
the plaquette operator [^]. Therefore, we pick b > 0.6. 

— b cannot be too big since we have to be below the deconfining transition for 
L = 5. Therefore, we pick 6 < 0.8 fll. 



• We use smeared space-like links and unsmeared time-like links. 

• We use the tadpole improved coupling bj = be{b) to set the scale and consider K xT 
Wilson loops W{K,T) with 1.5 < ^ < 12.5. This amounts to expectation values 
of Wilson loops that range from 0.82 to 2 • 10~^. 

• Keeping K fixed, we fit 

lnW{k, t) = -a- m{k)t; (2.1) 

where k = ^ and t = ^ are the dimensionless extent in the space and time direc- 
tion respectively. m{k) is the dimensionless ground state energy. This fit assumes 
that there is a perfect overlap with the ground state. Note that a should be zero 
since W{k,Q) = 1. Any small deviation from zero seen in the fit is due to the 
contribution from excited states. This can be seen by noting that W{k,t) should 
behave as e-'^e-""('=)* + (1 - e-'')e-™i(^)* with mi{k) > m{k). Then, lnW{k,t) = 
—a — m{k)t + ln\l + (e"" — 1) Q~["^i(k)-m{k)]f^ _ rpj^^ term is numerically insignifi- 
cant in the range of t being considered. 

• Finally, m(A;) is fit to ab'jk + cobj + The combination ^/abJ is plotted as a function 
of bj'^. We expect lattice spacing effects to lead off as bj^ in Yang- Mills theories and 
this is indeed the case in Fig. ||. The continuum limit extracted from this figure was 
quoted in Eqn.(|L3[). 



The use of smeared links improves the measurement of Wilson loops. They enhance 
the overlap of the space-like sides of the Wilson loops with the ground state. This increases 
the signal relative to the fiuctuations and simplifies the t behavior of the loops 0. One 
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step in the iteration takes one from a set uj,^\xi, X2,t) to a set U^^^'^' {xi, X2,t), by the 
fohowing equation: 

X['+'\xi,X2,t) 



(1 - f)Ui'\x) + ^U^\xi,X2,t)Ui'\xi,X2 + l,t) M\xi + l,X2,t) 



Xi'+^\x,,X2,t) 



t 



2 

?72^*^(xi,X2 - l,t) ' U^''{xi,X2 - l,t)U^''{xi + l,X2- l,t) 

(1 -/)[/« (x) 
2 



+ ^[/{'^(xi,X2,t)C/f (Xl + l,X2,t) \ui'\xi,X2 + l,t) 



t 



2 



?7f^(xi - l,X2,t) ' U^''{xi - l,X2,t)U'l''ixi -1,X2 + l,t) 



t 



r(0 



C/^'+')(xi,X2,t) =X<,'+^)(xi,X2,t) 



[4*+')(xi,X2,t)]tx(^+^)(xi,X2,t) 



; k = l,2 (2.2) 



Note that time-hke hnks, U3{xi,X2,t), are not smeared. Also note that smearing only 
involves space-like staples. There are two parameters, namely, the smearing factor / and 
the number of smearing steps n. Only the product r = fn matters, and / plays the role of 
a discrete smearing step. For a given r, the overlap of the smeared loop with the ground 
state does not depend on / as long as it is small. But the overlap of the smeared loop 
with the ground state does depend upon r. We set the value of the smearing parameter to 
r = 2.5 by choosing / = 0.1 and n = 25. To study the effect of varying r, we also consider 
r = 1.25 (/ = 0.05 and n = 25) at one coupling. 



3. Extraction of string tension 

SU{N) gauge fields were generated on a 5^ periodic lattice using the standard Wilson 
action. One gauge field update of the whole lattice Q is one Cabibbo-Marinari heat-bath 
update of the whole lattice followed by one SU{N) over-relaxation update of the whole 
lattice. A total of 1500 such updates were used to achieve thermalization. Measurements 
were separated by 10 such updates and all estimates are from a total of 832 such measure- 
ments. Errors in all quantities at a fixed b and N were obtained by jackknife with single 
elimination. 

The ground state energy m{k) obtained as a function of /c = ^ is fit to 

m{k) = abjk + cobi + ^ (3.1) 

k 

We expect ab'j to approach a finite value in the continuum limit {bj oo). The same is 
expected for ci. For large, unsmeared or symmetrically smeared Wilson loops with t ^ k, 
the universal value is — ^ ~ —0.13 ||8| rather than the value — ^ for Polyakov loops |9| that 
was seen in [||. The cobj term is present due to the perimeter divergent contribution, and 
therefore it is expected to logarithmically diverge in the continuum limit. Since, we do not 
smear the time-like gauge fields, the divergence in this term is not tamed. 

The method will encounter difficulties in extracting the physically relevant string ten- 
sion from Eq.( |3.1[ ) if cobj is large. However, because we do not go to very weak couplings. 
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we see in Fig. |2| that cqIj is not too large. It is not necessary to go to weaker couplings 
since the string tension computed at the couplings we have chosen can be used to get a 
good estimate of the string tension in the continuum limit as is evident in Fig. ||. 

The three parameter fit of m{k) as a function of k is shown in Fig. The fit has two 
degrees of freedom at the coarse lattice spacing of 6 = 0.6 and has four degrees of freedom 
at the fine lattice spacing of 5 = 0.8. As mentioned before, errors in ab^, ci, and cobj are 
obtained by jackknife with single elimination. Unlike the estimate of the leading coefficient 
abj, the estimates of the sub-leading ones are not as reliable. Consider the dotted line 
and dot-dashed line in Figure |2|. The constant term in each of these lines is obtained by 
evaluating cobj + ci/k at k = 12.5. The coefficient of the linear term in each is set to the 
same value as the one in the full three parameter fit. A comparison of the dotted line with 
the solid line and a comparison of the dot-dashed line with the dashed line shows that the 
1/k term becomes relevant only if A; < 7. There are only two data points at 6 = 0.6 and 
three data points at 6 = 0.8 where the 1/k effect is significant. As such, we do not expect 
our estimate of ci and cobj to be as reliable as the estimate of ab'j. 

The behaviors of cobj and ci as a function of bj^ are shown in Fig. ^ and Fig. |^. 
Between these two terms, cobj is the dominant one. The rise in cobj at smaller bj^ is 
consistent with the presence of a In 6/ term from a perimeter divergence. The estimate 
of ci is least reliable since it is dominated by the other two terms in the fit of m{k) as a 
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function of k. A fit of ci as shown in Fig. ^ is consistent with the existence of a continuum 
hmit. One should note that the number of degrees of freedom in the three parameter fit 
of m{k) increases as bj decreases and this will have an effect in the determination of the 
sub-leading terms. We believe this is reason the fit of ci versus bj does not pass through all 
the data points. Since we do space-like but not time- like smearing and since our loops do 
not generally have t ^ k, it not surprising to see a result that disagrees with the universal 
value but has the same sign and order of magnitude. 

4. Extraction of m{k) 

The dimensionless ground state energy m{k) is extracted at a fixed k by fitting \nW{k,t) 
to —a — m{k)t as discussed in Sec. While m{k) should be independent of the smearing 
parameter r = fn, the value of a is expected to depend r. 

We will use b = 0.8 as the coupling to illustrate the extraction of m{k). Figure |5| and 
Fig. ^ show the performance of the fit for two different values of r, namely, 2.5 and 1.25 
respectively. The solid circles show the data points without errors. The solid lines show 
the fit of the data. Seven values of t were used to fit the data at one k, and data at seven 
different values of k were fitted. This amounted to all Wilson loops from 1 x 1 to 7 x 7 
on the 5'^ lattice. The set of thermalized configurations used at r = 2.5 is statistically 
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-0.0'«i(6)-0.0184(15)/b' 



Figure 3: The behavior of the coefRcient 
cobi in the fit of m{k) vs k is consistent 
with the presence of a In 6/ term due to the 



Figure 4: The behavior of the coefRcient ci 
in the fit of m{k) vs k shows the existence of 
a continuum hmit. 




Figure 5: Plot of lnW(k,t) as a function 
of t for seven different values of fc at 6 = 0.8 
with r = 2.5. 




Figure 6: Plot of lnW{k,t) as a function 
of t for seven different values of fc at & = 0.8 
with T ~ 1.25. 



independent from the set used at r = 1.25. The fit parameters are shown in Table || and 
Table H Only the average values of the fit parameters are listed. 

Investigation of Table || and Table § shows that m(A;) does not depend on r. There is 
a small difference in the two values of m{k) at a fixed k for the two different values of r if 
k is large. But Fig. ^ shows that this difference is within errors. Furthermore, the fitted 
values of ab^ for the two different values of r are the same within errors. 

k 1.62 3.23 4.85 6.47 8.08 9.70 11.31 
a 0.001 0.003 0.009 0.019 0.055 0.047 0.071 
m{k) 0.133 0.218 0.286 0.347 0.399 0.464 0.517 

Table 1: Fit parameters corresponding to the fit lnW{k,t) — —a — m{k)t for seven different values 
of fc at & = 0.8 with r = 2.5. 

The values of a in Table [l| and Table |^ do show a variation with r and k. Since a 
smaller value of r implies less smearing, the overlap with the ground state is less for smaller 
r, and this results in a larger value of a at smaller r. The value of a is very close to zero 
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for small k indicating excellent overlap with the ground state for the chosen value of r. As 
k increases, the length of the loop increases and the perimeter divergence has a stronger 
effect. This results in a larger value of a as increases at a fixed r. 

k 1.62 3.23 4.85 6.47 8.08 9.70 11.31 
a 0.002 0.012 0.029 0.054 0.102 0.114 0.144 
m{k) 0.133 0.218 0.287 0.349 0.404 0.468 0.526 

Table 2: Fit parameters corresponding to the fit In W{k, t) = —a — m{k)t for seven different values 
of fc at 6 = 0.8 with r = 1.25. 




5. Finite N effects 



Two issues need to be addressed with the analysis performed so far. We have fixed our value 
of A'^ assuming finite effects are small. If A^ is not large enough, finite A effects need to 
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be addressed. In addition, we also have to address finite volume effects since continuum 
reduction is valid only in the N ^ oo limit. 

We expect m{k) to have a fixed limit as ^ oo at a fixed k, L, b and r. Indeed, this 
is the case as shown in Fig. ^ where the results for m(k) as a function of k are shown for 
b = 0.8 with T = 2.5 on 5^ lattice. All three fit parameters are consistent within errors all 
the way from = 23 to = 47. The only glitch one sees is at A; ~ 8. This corresponds to 
K = kbi = 5, which is the linear extent of the lattice. One can argue that there are larger 
finite N effects at strong coupling for K = L. Since the fit of m{k) involves several values 
of k, the larger effect at this particular value of k is diminished in the extraction of ab'\. 




Figure 8: The ground state energy m{k) as a function of k for five different values of TV at 5 = 0.8. 



Since finite effects can be ignored at A^ = 47, we also expect there to be no appre- 
ciable finite volume effects at this value of A^. This point is illustrated in Fig. |9| where the 
result for m{k) is plotted at & = 0.6 and r = 2.5 on 4'^ and 5'^ lattice. We used b = 0.6 for 
this comparison since we have to be in the confined phase on 4^ lattice. Figure ^ shows 
that the two values of m{k) at a fixed k are consistent with each other within errors. The 
same is the case for the fit parameter (t6|. This is not the case for ci and c^bj, and this is 
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probably due to a three parameter fit using only five data points. Sub-leading coefficients 
are expected to depend sensitively on the data points. Since we are primarily concerned 
with the value of the string tension in this paper and since all our results are based on data 
taken on 5^, we expect the final result to be free of finite and finite L errors. 
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Figure 9: The ground state energy m{k) as a function of k on two different lattices at b = O.i 



6. Creutz ratio 

It is natural to ask how the Creutz ratio ll 



W(K,J)WiK-l,J-l) 

' ^ W{K,J -1)W{K -1,J)' ^ ' 

performs as an observable from which to extract the string tension. If we were to use 
Creutz ratios, we would have smeared all links using all staples. But one can still ask how 
the Creutz ratio behaves with the asymmetrically smeared links. We show this for square 



loops (J = K) at 6 = 0.8 and r = 2.5 in Fig. IC. The solid lines show the estimate for the 



^/o'bi as obtained from the analysis in this paper. There is no evidence for a plateau in 
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the Creutz ratio in the range of k shown in Fig. |T^. It is possible the situation would be 
different if we had smeared all links. 
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Figure 10: Behavior of the Creutz ratio for square loops at 6 = 0.8. 



Each data point in Fig. 10 is obtained using only four different Wilson loops, i.e. four 
of the data points in Fig. ^. This is quite different from the analysis in this paper. Seven 
different Wilson loops in Fig. ^ are used to extract one m{k) point in Fig. ^, and the loops 
used for different k form independent sets. Then the m{k) are fit to determine the string 
tension. Both folded and unfolded loops contribute together. This is the main reason we 
succeeded in extracting the string tension using the range of Wilson loops considered here. 
To extract the string tension using Creutz ratios, larger loops and therefore larger statistics 
and possibly larger N would be needed. 



7. Conclusions 



We used Wilson loops with smeared space-like links and unsmeared time-like links to 
obtain an estimate for the string tension in the large limit of three dimensional Yang- 
Mills theory. Invoking large N continuum reduction, we included Wilson loops larger 
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than the size of the lattice. Since we used smeared space-hke hnks, the Wilson loops 
for fixed length in space and varying length in time showed excellent agreement with a 
single exponential. The ground state energy so obtained was fit using three parameters to 
get an estimate for the string tension. The ground state energy exhibited short distance 
behavior at the shortest length used in the paper but large enough distances were used 
to get an estimate for the dimensionless string tension with small errors (Equation pT^ ). 
These results validate the method of continuum reduction for calculating quantities based 
on the space-time dependence Wilson loops. 
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